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Abstract 



We provide a fast and precise Mellin-space implementation of the 0{as) heavy flavor Wilson 
coefficients for charged current deep inelastic scattering processes. They are of importance 
for the extraction of the strange quark distribution in neutrino-nucleon scattering and the 
QCD analyses of the HERA charged current data. Errors in the literature are corrected. 
We also discuss a series of more general parton parameterizations in Mellin space. 



1 Introduction 



The scale evolution in Quantum Chromodynamics (QCD) can be performed very efficiently in 
Mellin space. There the evolution equations of the twist-2 parton distributions and the corre- 
sponding observables become ordinary differential equations, unlike in momentum-fraction space, 
where integro-differential equations have to be solved numerically. In Mellin space, moreover, 
the evolution equations can be solved analytically, provided that the anomalous dimensions and 
Wilson coefficients can be represented at complex values of the Mellin variable N, cf. [1]. The 
transformation mediating between both representations is 

M[f{x)]{N)= [\x x""-' fix) , (1) 

where f{x) may be distribution valued. The anomalous dimensions and massless Wilson co- 
efficients to 3-loop order are functions of the nested harmonic sums [2,3], see Refs. [4]. The 
nested harmonic sums Ss{N) are meromorphic functions in the complex plane with poles at 
the non positive integers and in case that the ffist k indices aj are equal to one, they diverge 
oc ln*^(A^) for large values of iV, [5]. They obey recursion relations for arguments N — )■ {N + 1) 
in terms of harmonic sums of lower weight, which allow for shifts parallel to the real axis. More- 
over, analytic asymptotic representations have been derived, [5,6], which are valid in the region 
of large values of A^, |arg(A^)| < n. One may derive effective numerical representations using 
the MINIMAX-method [7], cf. [8,9]. In this way the QCD observables at the level of twist-2 are 
represented analytically at general scales Q^, parameterizing the parton density functions by 
appropriate non-perturbative distributions at a starting scale Q^. This also applies for a wide 
class of other processes, cf. [10]. 

In many precision analyses heavy quark contributions have to be considered, which are more 
difficult to represent in Mellin space. For deep-inelastic scattering the contributions up to the 
(^(q;^) neutral current heavy flavor Wilson coefficients [11], which are available in semi-analytic 
form in momentum fraction space, were given in Mellin space in Ref. [12]. The running mass 
effects have been described in [13]. 

In the present note a fast and precise implementation of the 0{as) charged current heavy- 
flavor Wilson coefficients for deep-inelastic scattering is presented. In Section 2 we derive the 
corresponding Mellin-space representations for the charged current structure functions up to 
0{as)- As a significant part of the data emerges for large values of Q^, we also derive the repre- 
sentation in Section3 which is valid for ^ m^, with = —q"^, q being the four- momentum 
transfer, and the charm quark mass. We illustrate the precision of the Mellin-space imple- 
mentation by comparing to the structure functions in momentum fraction space in Section 4 . 
Finally we remark on the Mellin space implementations of various x-space distributions used in 
different analyses in the literature in Section 5, which may be important to account for a higher 
flexibility in the choice of the non-perturbative input distributions. This will allow to perform 
analogous analyses in Mellin space in the future. 



2 The Scattering Cross Section 

The charged current deep-inelastic scattering cross section for z/(z7)p — /(I) -|- X, l(J)p 
viv) + X in case of charm-quark production is given by 

^ ^ ^ y^2x¥l{x, Q\ m^) + 2 1 - y - y"—^ F^(x, , ml) 



dxdy 471 {M^ + Q^y\" IK ,^ , cj ' y " " Q 
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±[l-{l-yy]x¥l{x,Q',ml)j. (2) 

Here x = Q'^/yS,y = q.P/l.P are the Bjorken variables, / and P are the lepton and nucleon 
momentum, S = 21. P, Gp is the Fermi constant, Mw the mass of the W^-boson, M the nucleon 
mass, and Fj(a;,(5^) are the structure functions. The ± signs in (2) refer to incoming neutrinos 
(anti-neutrinos) or charged leptons (anti-leptons), respectively. Note that the scattering cross 
section does depend on the masses of the initial state particles, cf. [14]. Due to the inclusive 
kinematics the dependence on rric is only implicit. 

In the twist-2 approximation, invoking the parton model, one may decompose the nucleon 
wave function into individual partons and consider the excitation of single charm quarks in the 
transitions 

s' = s \V,,\^ + d \VJ^ ^ c, s' = s \V,,\^ + d \VJ^ ^ c , (3) 

with Vij the CKM matrix elements [15], and s and d the strange and down quark parton distri- 
butions. Here d and s denote the corresponding anti-quark distributions. For this transition the 
Bjorken variable x and the momentum fraction of the struck parton ^ are related by 

X = < A . (4) 

with A = Q'^/iQ^ + fn^). This phenomenon is called slow rescaling [16]. The heavy quark 
structure functions and parton densities can be expressed as a function of ^ G [0, 1], preserving 
Mellin symmetry. 

The 0{as) charged current heavy-flavor Wilson coefficients were calculated in [17] and cor- 
rected in [18] later. We will follow Ref. [18] and work in the fixed fiavor number scheme (FFNS). ^ 
We define 

JT = F5, T^ = ¥y2^, J'^,=n/^ (5) 
To 0{as) one obtains, after a Mellin transformation (1) over ^, 



(6) 



for exchange. In case of W~ exchange s' replaces s'. Here = a^/ (47r) denotes the strong 
coupling constant and g the gluon distribution. The massive Wilson coefficients for charged 
current deep-inelastic scattering in Mellin space are given by^ 



^The use of a variable fiavor number scheme needs care, since the scales at which one massive flavor can be 
dealt with as massless is process dependent. The corresponding scales are usually not fi^ vn? , cf. Ref. [19]. 

^Throughout the present paper we consider the scale derivative in the renormalization group operator as 
dld\xi{^?). 
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with 



(A,iV), z = l,2,3, 



(9) 
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g2 + ,„2^ ^ 



^S'M^^^^^l = 7^i?(^^)ln(^^-V^)+7P(A,iV) + /^l2(A,iV) (10) 



Here the leading order sphtting functions are 



p(°)(iV) = 8Tf 



2 A^(A^ + 1) 
N^ + N + 2 
N{N + 1){N + 2) 



25i(iV) 



/if (A, AT). (11) 

(12) 
(13) 



with Tp = 1/2, Cf = (iVc - l)/(2iVc) for ^f/(iVe) and N, = 3 for QCD. The function P(A, iV) is 
given by 



P(A,iV) 



M 

8Tf 



Pi?(z)ln 



1 - Xz 



(N) 



(1-A)z_ 

[A^ (iV2 + + 2) - 2A (A^2 + 2N) + 2N^ + 2N] 



XN{N + 1){N + 2) 
A {N^ + 3N + 2)- 2N^ 



P2(A,iV)--ln(l-A) 



XN^{N + 1){N + 2) 



(14) 



In the hmit A — 1 it takes the form 

r N'^ + N + 2 
P(1,N) = 8Tp <^ — -— ^^7— -- 
^ ' ^ \A^(iV + l)(Ar + 2) 

The functions hf^\X,N) read : 
/i?(A,A^) = 



In 



Si{N) 



N'^{N +1){N + 2) 



(15) 



4P,(A, iV) - [A^ (iV^ + iV + 2) - A (iV3 + 3N^ + AN) + 2N- + 2N] , ^ 



[X{9N^ + 8iV2 _ 5iV - 2) - 2iV2] 



AiV(A^ + 1) 

(1 - A) [A {N^ + 3iV + 4) - 2iV - 2] 
A2(iV + l) 



/i^(A,iV) = Cf 



4Pi(A,iV) + 



AiV2(7v + 1) 

[4A2 (iV2 + N)-X{N^ + 2iV2 + + 2) + iV^ _ ;y2] 



ln(l - A) 
(16) 



4 (iV2 + iV - 1) 



5i(iV)+45i2(iV) 



iV(iV + 1) 
(1 - A) [4X{N + 1) - + N] 



F2{X,N) 



N{N+1) 
[AX{N^ + N)- 9N^ -6N^ + N + 2] 



A(Ar+l) 



ln(l - A) 



(17) 



hl{X,N) 



.2..n (l-A)(iV2 + iV + 2) 



+AStiN) - 



N{N + 1) 

ln(l - A) - 



iV(iV+ 1) 
{3N + 2) (3Ar2 _ 1) 



/i?(A,iV) 



-4Tp 



A(iV + l) ' ' iV2(^iV+i) 

[A2(2iV3 + 3iV2 + 7N + 2)- X{2N^ + AN^ + 8N) + 2N^ + 2N] 



{N'^ + N + 2) 



-Si{N) 



2XN{N + 1){N + 2) 
(1-A) [X{N^ + N + 3) - N -1] 



(1^ 
F2{X,N) 



N{N + 1){N + 2) ' ' A2(A^ + l)(A^ + 2) 

A(2iV3 + iV2 + iV - 2) - 2iV2 



ln(l - A) 



/if (A, AT) 



-AT, 



2XN'^{N + 1){N + 2) 

[A3(2iV3 - QN^ + AN) + A2(-2A^3 + 7^2 _ ^ 3) - A(2A^2 ^ 4^)] 

2AiV(A^ + l)(iV + 2) 

' ^ A2(A^ + l)(iV + 2) ^ ^ 

A2(2Ar3 - 6Ar2 ^ _ ^(jv2 + SAT) - 2A - 2N^ 



(19) 



A(A^ + 2) 
(iV2 + jv + 2) 



-Si{N) 



hliX,N) 



N{N + 1){N + 2) ' ' 2XN'^{N + 1){N + 2) 

(A2(iV2 _ ;y ^_ 2) + A(2iV2 ^ 4JY) _ 2iv2 - 2iV) 



(20) 



-4T, 



2AA^(iV + l)(iV + 2) 



-F2(A,iV) 



A^(Ar+ l)(7V + 2) ^ ^ ^ A2(Ar + l)(Ar + 2) 
A(A^2^3iV + 2) -2iV2 



2XN'^{N + 1){N + 2) 
Here Sk{N) are the single harmonic sums [2,3] 



(21) 



(.-1)1 

and Fi 2(A, A^) denote the functions 

Fi{X,N) 
F2iX,N) 

r(i)-9 



^(fc-i)(]V +!) + (;,, A: eN\{0} 



(22) 



dx 

2; - 1 



ln(l - Xx) 



x^ -1 



dx 



1 - Xx 



(23) 
(24) 



We refer to H2 for Q = /i also in the form 



Hi^^'%X,N) = ATf< 



A - 2N{N - 3) - N{N^ + N + 2){2Si{N) + ln[A(l - A)]} 
2iV2(iV + l)(iV + 2) 



18(1- A) + 12(1- A)2 (l-A)2Fi(l,Ar,iV+l;A) 1 



(A^ + l)(iV + 2) 
+6A(1 - A) 



N 



N 



2Fi(l, TV + 1, AT + 2^ _ i(l, + 1, iV + 2; A) - 1 



(iV + 1) 



{N + l){N + 2) 



(25) 



to correct Eq. (70) in Ref. [20]. Here, 2-^1 denotes the GauB' hypergeometric function 

,F^{rJ, N + a + l,N + a + P + 2;0= dzz'^+^il - z)^{l - Cz)-\ 

B{N + a + l,/3 + 1) Jo 



(26) 



k\X 


0.1 


0.5 


0.9 





_l_q qqqqqqqqqqqqqqqqqq 




qqqqqqqqqqqqq5454T 


+ 1 111110872017648708 


1 


+0 821423460E-22 




1 9751 23n84486fi87F,-1 f) 


+0 000098847800695649 


2 


-0 050000000000000000 

\J • V./ \J \J \J \J \J \J \J \J \J \J \J \J \J \J \J 


-0 


250000001423445310 


-0 456776649795426910 


3 


-0.001666666666666666 


-0 


041666626227537189 


+0.048001787041481484 


4 


-0.000083333333333336 


-0 


010417267014473893 


-2.639926313632337045 


5 


-0.499999999997477753E-05 


-0 


003119649969153064 


+21.74650303005800508 


6 


-0.333333333484935779E-06 


-0 


001072643292472726 


-119.2024353694977123 


7 


-0.238095231857886535E-07 


-0 


000249816388951957 


+442.8630812379568060 


8 


-0.178571609458723466E-08 


-0 


000477192766818736 


-1147.492168730891145 


9 


-0.138885135219891312E-09 


+0 


000607438245655042 


+2095.900562359927839 


10 


-0.111166985797600287E-10 


-0 


000939620170907808 


-2688.830619105423874 


11 


-0.903191209860331475E-12 


+0 


000876266420390668 


+2371.811788035594662 


12 


-0.800460593905461726E-13 


-0 


000571683728308967 


-1370.383159208738358 


13 


-0.439528295242656941E-14 


+0 


000219008887781703 


+467.1985539745224063 


14 


-0.107994848505839019E-14 


-0 


000041032031472769 


-71.30766062862566906 



Table 1: The expansion coefficients a^, Eq. (33) for special values of A. 

Many of the contributing functions have known Melhn transforms given before in [2,3]. Their 
analytic continuations to complex values of iV are known, cf. [5,6,8,9]. The new functions (23,24) 
are solely related to integrals of the kind (26) for G N. They obey the following recursion 
relations 



Fi(A,iV + l) 
F2(A,iV + l) 



i^i(A,iV) + 



A 



A^+ 1 



1 
A 



F2(A,iV) 



F2(A,Ar + l) 
1 1 - A 



iV+ 1 



ln(l - A) 



(27) 
(28) 
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The singularity structure of Fi 2(A, A^) for iV G C can be seen using the representation 

°° A' 

F,{\,N) = - J2 [^(l + N + l)-ij{l + l)] (29) 
1=1 ' 

°° A' 

F2(A,iV) = -A^Vt : --, r, (30) 

^ ^ ^ (iV + / + !)(/ + !) ' ^ ' 

where %jj^''\z),k > denotes the polygamma function. Both functions possess poles in at 
negative integers. 

For -Fi(A, N) one may derive a sufficiently precise representation using the MINIMAX-method. 
The function Fi can be written as 



(31) 



Fi{\,N) = -!^^-N£{\,l) + J2ai {N-l-l[S,{N + l-l)-S,m'^ , 

S{\,x) = ^(l-Ax)(l-ln(l-Ax)) . (32) 
A 

The function S under the integral is approximated by the adaptive polynomial, cf. Table 1, 
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S{X,x) = Y,ak{X) x'' . (33) 



fc=0 



Knowing the difference equations (27,28) one may shift Fi 2 parallel to the real axis. Usually one 
attempts to shift towards the asymptotic region and applies an analytic representation there. 
However, as sometimes recommended, one may also turn the view [21], and rather consider 
-^1,2 (A, A^) inside the unit circle, to obtain an even more beautiful result : 

F,{X,N) = -^[Sl{N) + S,{N)] 
00 

- {LW2(A) - + i^i,fc+i(A) - Hi,k+i{l)} (34) 



k=l 



1 f ~ AT 

^2(A,iV) = T E(-^)'[LWi(A)-l]-^^ 
I k=i 



(35) 



for |A^| < 1. Here, Lifc(a;) denotes the classical polylogarithm [22], Hg{x) a harmonic polyloga- 
rithm [23], 

H,,,{x) := r dz ^^M^ , (36) 
Jo V — z 

and Cfc Riemann's (^-function at an integer argument. In the limit A ^> 1 one obtains 

Fi(l,Ar) = -1\S\[N) ^ S^{N)\ (37) 
F2(l,Ar) = -S^{N) . (38) 

The serial representations (34, 35) are quickly converging and deliver even more pre- 
cise results than using (31, 33). The contour integral is performed along the line 
(— 00, — a), (c, — a), (c, a), (— 00, a), with a = 0.5. The choice of the parameter c depends also 
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on the rightmost singularity of the non-perturbative distribution f{x) and was chosen with 
c = 1.5 in the figures given below. The inverse Mellin transform is obtained by 

f{x) = ^!^-J \m [x-'-'^Mlfjis + ia)] ds + \m [ix-''-'M[f]{is + c)] ds^ . (39) 

We extended the integral to 1000 units parallel to the real axis. The recursion relations (27, 
28) together with the series (34, 35) allow to compute -Fi 2(A, A^) along the integration contour. 
Here usually the first 30 terms in the infinite sums (34, 35) provide sufficient accuracy. 



3 The Massive Wilson Coefficients in the Asymptotic Re- 
gion 

In many applications the value of A becomes close to one. Within this kinematic region Q"^ ^ 
the Wilson coefficients take a simpler form, and the functions needed to express them are only 
harmonic sums. This has been shown up to 3-loop order in case of neutral current deep-inelastic 
scattering in [24]. At 0{as) one obtains the following relations : 



Here, A^^^ denote the 1-loop massive operator matrix elements derived in the neutral current 
case and C'fg^g) are the massless 1-loop Wilson coefficients, [25]. Eqs. (40, 41) are derived by 
expanding (16-21) for A — )■ 1. At 0{as) there is no pure-singlet contribution. It turns out, that 
A^q^'^^'^ vanishes as expected, because closed massive fermion loop contributions can occur at 
O(a^) earliest. In the charged current case the interaction transmutes the massless s'-quark into 
the massive c-quark. The massless quark-loop contributions to Aq^, occurs with a combinatorial 
factor 1/2 compared to the neutral current case. These terms vanish, because the corresponding 
diagrams are scaleless. Because of this no quarkonic operator matrix element contributes to (40). 
Eqs. (40, 41) agree with Ref. [26]. 

The massless Wilson coefficients obey 

C., (iV, ^) = 5,, + f: o^Cj^ {n, ^) , (43) 



where 



^'■'(^.f) = (44) 

^■'(^•f) = 5n'?(iV)lng)+cS(iV) (46) 
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lO"* 10"^ 10"^ 10"'' 10° 

X 

Figure 1: Relative accuracy of the charged current structure functions Fj to 0{as) at = 
10,100,500 GeV^, comparing the implementation in Mellin- and z-space. Here both for Fi 2(A,A^) 
the analytic representation (34, 35) and their recurrences were used. 



and [33] 



.(1) 
-3,q 



(N) 



(iV) 



.(1) 

-2,!? 



{N) 



-2,g 



Cf 

2,g 
''2,3 

Tf 
. 



(3iV2 



iV+ 1 
3A^-2) 



(N) - C 



N{N + 1) 

2(2iV + 1) 



SiiN) + 2S^,iN) - 2S2iN) - 



{9N^ + 2N^ -5N-2) 
iV2(iV+ 1) 



(N) 



N{N + 1) 

^iN + l){N + 2) 
A{N^ + N + 2) 



N{N + 1){N + 2) 



Si{N) 



4(iV3 _ 4JY2 _ _ 2) 
m{N + 1){N + 2) 



(46) 

(47) 
(48) 

(49) 

(50) 
(51) 



We note that the global sign for Eq. (A1.17) of [26] has to be reversed. The gluonic contribution 
in charged current heavy quark production has been calculated before using two finite quark 
masses for e~p scattering in Refs. [27-29] and both for e~N and e~^N scattering in [29]. For 
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calculations in case of neutrino-nucleon scattering see [30]. We repeated the calculation of the 
gluonic contribution for a massless s-quark in the MS-scheme. 
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Figure 2: The charged current structure functions for charm production for different values of Q'^ , 
using the ABKM09 parton parameterizations [31]. We indicated the kinematic range being accessible 
at HERA, (a)-(c) : the structure functions for exchange, (d)-(f) difference between the structure 
functions for and W~ exchange. 



4 Numerical Results 

We compare the Mellin space representation given in Section 2 with the representation in x-space 
of Ref. [18] using the reference distribution 

x]{x) = x-^-^{\-xf (52) 
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0.02 
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X 

Figure 3: The difference of the structure function F2 for exchange for an isoscalar and a proton 
target using the ABKM09 parton distribution functions [31] in depending of a; and Q"^. 

for both the quark and gluon densities and determine the relative accuracies |(5Fj|/Fj for different 
values of in the massive Wilson coefficients choosing rric = 1.5 GeV and the corresponding 
values of aNLO(Q2) p^]^ 

If one employs the MINIMAX representation relative accuracies of better than 3 x 10~^ are 
reached below x = 0.5. For x > 0.5 the relative accuracy becomes worse. In this region, 
however, the charm contribution is very much suppressed, as shown in Figures 2 and 3. Using 
the representations (34, 35) the relative accuracy is improved and amounts to 5x 10~^ at x = 10~^ 
or better growing to ~ 10~^ for x ~ 0.4, cf. Figure l,a-c. Beyond this value the relative accuracy 
becomes worse, but also the sea quark distributions are very small in this region. For comparison 
we note that in [20] accuracies of 0.015 to 0.002 were obtained. 

In Figure 2 we show the charged current heavy flavor structure functions 2xFi, Fi and XF3 for 
ep-scattering and the kinematics at HERA for exchange and the difference of the structure 
functions for and W~ scattering using rric = 1.5 GeV, referring to the ABKM09 NLO 
distributions in the fixed flavor scheme [31]. Results are shown for = 10, 100,500 and 1000 
GeV^, (a^LO = 0.2399,0.1666,0.1379,0.1283). The dominant contributions are those due to 
ly^-gluon fusion, which is the same for and W~ exchange. All structure functions rise 
towards small values oix. The difference between the and H^~-exchange structure functions, 
on the other hand, receives its main contributions in the valence region with smaller scaling 
violations as in case of VT^-gluon fusion. 

In Figure 3 the difference of the structure function F^'^ for an isoscalar and proton target is 
shown. This also is a valence-like function. With the present precise Mellin-space implementa- 
tions at hand one may perform QCD fits including the charged current heavy flavor contributions 
in an efficient way. 

5 Parameterizations in Mellin Space 

Parameterization of the non-perturbative input distribution of the parton distribution functions 
at a starting scale /Iq containing sufficient flexibility is an important condition for the QCD 
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analyses of deep-inelastic and hard scattering data. In the literature different parameterizations 
of distributions are used. A wide-spread shape has the form 



K 



Mx,fil) = -x)^P(x), P{x) = Y.lkx''^ . (53) 



k=l 

The corresponding Mellin transform reads 

K 



/i(iV,/i2) = A^7fc5(iV + a + 5fc,/3 + l) , (54) 



k=l 



with parameters such that N + a + 5^, P + ^ 7^ — /, V/ G N and B{x, y) = T{x)T{y)/T{x + y). 
Eq. (54) is analytic under this condition. 

One may adopt the attitude to represent the non-perturbative momentum distributions of the 
partons, which are L'^[0, 1] measurable functions, in terms of orthogonal polynomials, see [32,33]. 
If the argument of these polynomials is x or a real power of x one may refer to structures as (54), 
provided the Mellin transforms exist. The fastest convergence is obtained choosing Laguerre 
polynomials [34] of the argument ln(l/x), [33]. 

f i{x, fil) = ^CkLk (\n(-jj , (55) 

with 



^.W = E(/)(-i)'7r- (56) 



The Mellin transform reads 

°° /AT 1 \ fc 



N ^ ^ \ N 

k=0 ^ 

However, also other parameterizations are used in x-space analyses. One of them [31,35] 
refers to the following shapes 

/,(x,/x2) = Ax"(l-x)V^(i+''^), (58) 

partly with 6 = 0. Let us first consider the simplified case 6 = /3 = as an illustration. Here, 
one obtains 

which is easily recognized as a quickly convergent series in k. In the general case we obtain 

^{N + a + k + l,/3 + l). (60) 



fc=0 1=0 ^ ' 



The derivatives of the Beta-function in (60) can be calculated in a recursive way, see [3]. Let 

iJ^i?(iV + a, /3) = i?(iV + a, /3)Afe, Aq = 1 . (61) 
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Then one obtains 

1 ^ 

Afc = -^5^Afc_^/, , (62) 
1=1 

with 

= (TTI)! [^^'-'\N + a) - ^p^'-'\N + a + (3)] . (63) 

Yet another parameterization was used in [36], 

f\{x, = Ax^'il - x)^e^%l + 6xy. (64) 

We first consider the case 5 = 0. The MeUin transform of (64) is then given by 

fi{N,Q^) = AB{N + a,P + l) iFi{a + N,f3 + a + N + l--f) 

= A 5(iV + a,/3 + l) e%Fi(/3 + l,iV + a + /3 + l;-7) , (65) 

where iFi{a, P; z) denotes the confluent hypergeometric function [37]. It obeys the recursion 
relation [38] 

iFi(a, P + l;z) = - . [/3(1 - /3 - z),F,{a, /3; z) + /3(/3 - l)iFi(«, /3 - 1; z)] . (66) 

2;(p — a) 

For large values of the asymptotic relation [38, 39] 

,F,(/5 + i7V + a + /3 + i--7) = y r(/^ + i + ^)r(iv + « + /^ + i)(-7r 

+0(|iV + a + /3 + l|-") (67) 

holds. Eqs. (66,67) may thus be used to compute the analytic continuation of (65). 
For the additional factor 

(1 + fe)" (68) 

one may extract an integer power from 77 such that the remainder,-?/, obeys f] G [—1/2, 1/2]. The 
integer contribution fj = rj — f] is modifying the above representations but leads to the same 
structure. The remaining factor possesses the well converging representation 

(l + 5x)^ = f;(^)5'a;', (69) 
and leads to a corresponding generahzation of relations (65, 67). 



6 Conclusions 

A fast and precise Melhn-space implementation of the 0{as) massive charged current Wilson 
coefficients for deep-inelastic scattering is provided. This process is of relevance for the un- 
derstanding of di-muon production in deep-inelastic neutrino-nucleon scattering and for the 
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deep-inelastic charged current data measured at HERA in particular w.r.t. the extraction of 
the unpolarized strange-quark distribution. To represent the Wilson coefficients we used both 
the MINIMAX-method and completely analytic representations and exploited recurrences in the 
complex- valued Mellin variable for N ^ N + 1. In wide kinematic ranges relative accuracies 
of 10~^ and better are obtained for the deep-inelastic structure functions. Errors in the litera- 
ture have been corrected. We present phenomenological applications for the structure functions 
Fj(a;,Q^), i = 1,2,3 both for and W~ exchange. Furthermore, we provided Mellin-space 
representations for a wide range of representations used in different x-space analyses. This allows 
both for more flexible choices of the parton distributions at the input scale using Mellin-based 
codes in data analyses and for comparisons with fits given in the literature. The corresponding 
code is based on ANCONT [8] and available on request from Johannes.BluemleinOdesy.de. 
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